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1 Introduction 



The computation of the effect of a perturbation in the potential V(x) of a one- 
dimensional Schrodinger operator 

H = -d 2 x + V(x) , (1.1) 

defined on the whole real line on the continuous part of the spectrum of H is a 
standard (albeit non trivial) procedure in perturbation theory. For example, the 
response of the reflection amplitude r(k) at momentum k (in a setting where there is 
source at x — +00) to an infinitesimal change in the potential is 

8r(k) _ (t(k)<Kx,k)) 2 

SV{x) 2%k 1 ' ' 

The derivation of this result is given in the Appendix. Here t(k) is the transmission 
amplitude and t(k)(j)(x, k) is the solution of the Schrodinger equation 

(-d 2 x + V{x) - k 2 ) t>{x, k) = (1.3) 

which satisfies the scattering boundary conditions of this problem: 

t(k)(f)(x,k) = t(k)e~ lkx + o(l) x-^-00 

t(k)(f)(x,k) = e' lkx + r{k)e lkx + o(l) x -> 00 . (1.4) 

Obtaining the kernel inverse to (1.2), i.e., the response of the Schrodinger potential 
V(x) (as well as the response of the wave function) to a change in the reflection 
amplitude r(k), is a much more difficult problem, which we solve in this paper. 

The explicit formulas we derive in this work express the local response of the 
potential and wavefunctions (see (3.13) and (3.15)) to a change in the reflection 
amplitude. Thus, our results add to the information which can be gleaned from the 
well-known trace identities of [1], the lowest of which reads 



oc 

/ 



V(x) dx = — I log [l - |r(A;)| 2 l dk - 4 f>* , (1.5) 

00 



00 N 
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where E\ = —nf are the iV bound state energies (with ki > 0), which tell us only 
about that response integrated over space. 

A possibly interesting application of the results of this paper might be the inves- 
tigation of perturbations around reflectionless potentials with arbitrary numbers of 
bound states[2], which play an important role in supersymmetric quantum mechanics[3] 
and the theory of solitons. 

It would be useful at this point to introduce some additional notations and recall 
some basic facts, which will be used later on. For k real, H — k 2 is real, and therefore 
W*(x, k) = ty(x, -k) for any solution of (1.3). It follows that 

r*(k)=r(-k) and t*(k) = t(-k) (1.6) 

in (1.4). 

^(x, k) and ^*(x, k) are linearly independent solutions of (1.3), and the continuous 
spectrum is doubly degenerate at each k 2 > 0. In particular, <f>(x, k) and (f)*(x, k) form 
a basis. Since they are degenerate in energy, their Wronskian is a non-zero constant. 
Equating its values at x — > ±oo, we obtain the probability conservation relation 

|r(A0| 2 + |f(fc)| 2 = l. (1.7) 

An equally suitable basis is the pair of solutions ip(x, k) and ift*(x, k) of (1.3), in which 
if){x, k) obeys the boundary condition 

ifj(x,k) = e- ikx + o(l), x-^oo. (1.8) 

Thus, ip*(x, k) corresponds to a setting in which there is a source at x — — oo. 
We see from (1.4) and (1.8) that 

(p(x,k)- 



1 : ^(x,k) + ^-i;*(x,k) 



> t(k)^ x ' ' t(k) 

This relation must hold for all x, since ip(x, k) and ip*(x, k) form a basis everywhere. 
Adding to it the linear combination for <p*(x, k), we may write the relation between 
the two bases as 



/ (j)(x, k) 
V <t>*{x,k) 



( 1 r(k) \ 

' t(k) t(k) 

r*(k) 1 

\ t*(k) t*(k) I 



( ip(x, k) \ 
\ ^*{x,k) J 



1.9) 



Note that the transformation matrix has a unit determinant. The inverse transfor- 
mation is thus 



/ ip(x, k) 
V i>*(x,k) 



I _J_ r ( k ) \ 

I t*(k) t(k) 

r*(k) 1 

V t*{k) t(k) J 



4>(x,k) \ 
<j>*{x,k) j 



(1.10) 



This paper is organized as follows: In the next section we present a lightning review 
of inverse scattering theory. In particular, we discuss the Gelfand-Levitan-Marchenko 
equation and its properties. We show that its solution is simply the boundary column 
of its resolvent kernel. 

In section 3 we compute the variational derivative of the solution of the Gelfand- 
Levitan-Marchenko equation with respect to the reflection amplitude. Then we derive 
from it the corresponding derivatives of the potential and wavefunctions (Eqs. (3.13) 
and (3.15), respectively) in closed form. 

In section 4 we demonstrate the consistency of our results by comparing their 
integrated form against known facts. 

Finally, in the Appendix we provide some useful technical details. In particular, 
we present the derivation of (1.2), and also discuss briefly the case of reflectionless 
potentials. 

2 The Gelfand-Levitan-Marchenko Equation and 
Its Solution 

According to inverse scattering theory (1ST) [4, 5, 6, 7], a Schrodinger operator (1.1) 

oo 

whose potential V(x) tends asymptotically to zero fast enough, such that J \V(x) | (1+ 

— oo 

I a; I) dx < oo, and thus supports only a finite number N of bound states, is uniquely 
determined by the so-called scattering data. The scattering data are the reflection 
amplitude r(k), and a finite set of 2N real numbers 

Ki > k 2 > ■ ■ ■ k n > and ci, c 2 , • • • cn , (2.1) 

where E\ = —nf is the Zth bound state energy, and where q appears in the asymptotic 
behavior of the /th normalized bound state wave function as Qexp — Ki\x\ (and thus 



determines its "center of gravity"). 

The prescribed reflection amplitude r(k) can be taken as any complex- valued 
function which satisfies (for k real): 

r(-Jfc) = r*(k) 

\r(k)\ < 1, A;^0 and 

r (fc) = O(j) , |fc|-oo. (2.2) 



The first two conditions were already mentioned in (1.6) and (1.7), and the third one 
reflects the fact that V(x) is a small perturbation at high energy. In addition to these 
conditions, there is a less obvious technical condition that the Fourier transform 
B(x) = J r Me ikx dk, should satisfy the bound 7(1 + \A)\^\dx < oo. Note 

— oo — oo 

further that t(k) is completely determined by r(k) as 



- |r(£;)| 2 \fJ{k-iKi) \2ni q-k 

Given the scattering data, 1ST instructs us to determine a certain real transfor- 
mation kernel K(x, y), bounded on the domain y > x, which maps the wave functions 
of the free Schrodinger operator H = —d 2 x onto those of the operator H in (1.1). For 
example, the left moving wave e~ tkx is mapped onto 

oo 

tj}{x, k) = e~ lkx + J K(x, y) e~ iky dy , (2.4) 

X 

which is evidently the solution of (1.3) satisfying the boundary condition (1.8) men- 
tioned above. Finally, the potential V(x) in (1.1) is determined by K(x,y) according 
to 

V(x) = -2-£-K(x,x). (2.5) 

The kernel K(x, y) is determined as the solution of the Gelfand-Levitan-Marchenko 
equation (GLM) [4, 5, 6, 7], 

oo 

K(x, y) + F(x + y) + J K(x, z)F(z + y) dz = , (2.6) 



where the real function F(x) is 



N 1 7 

= X! c « e""'" + ^ J r(k)e lkx dk . (2.7) 

For fixed (2.6) is a Fredholm integral equation of the second type in the unknown 
function 3>(y) = K(x,y), 

oo 

*(v) = f(v) + ><fN(y,z)*(z), (2.8) 

with symmetric real kernel and given function 

7V(y, = N(z, y) = -F(y + z) , f(y) = -F(x + y) = N(y, x) (2.9) 

respectively, and spectral parameter A = 1. 

It is known that (2.6) has a unique solution, i.e., the Fredholm determinant of 
(2.8) is not null at A = 1 [5, 6]. It is easy to demonstrate this property in the case of 
reflectionless potentials (for which r(k) = for all k), as we show in the Appendix. 

The unique solution of (2.8) at A = 1 is given by 

oo 

*(v) = /(v) + / R(y,z;l)f(z)dz, (2.10) 

X 

where R(y,z;X) is the resolvent kernel of (2.8). Note that for real values of A, 
R(y, z; A) is manifestly real, when it exists. 

It is useful at this point to introduce the operator iV and the vectors |$) and 
|/), which correspond to the kernel N(y,z) and functions and f(y). Thus, in 
obvious notations, N(y,z) = (y\N\z) = —F(y + z), (y\f) = f(y) = (y\N\x) and 
®(y) = (y\®)- Then, it is easy to see from (2.8) that 

|$) = — UH/> = — \X). (2.11) 

1 — N 1 — N K J 

Similarly, from (2.10) we deduce that R(y,z; 1) = (y\R\z), where 

1 N 
R=— -^-1 = —^. (2.12) 
1-N 1-N 



oo 



Thus, by comparing (2.11) and (2.12), we conclude that |<3>) = R\x) , i.e., 

*(y)=K(x,y) = (y\R\x). (2.13) 

The solution of the GLM equation (2.6) coincides with the xth column of its resolvent 
kernel. It is manifestly a real function of x and y. 

3 The Variational Derivatives 

In view of (2.13), it is straightforward to compute the variation of K(x, y) under small 
perturbations in N. Thus, consider a perturbation N — > N + 8N, which induces the 
variation R ^R+^SNj^. Consequently s -§*$- = + £)| a )<&|(l + R)\z). 
Thus, from (2.13) 

= (5(y - a) + R(y, a; 1)) (S(b - x) + K(x, &)) . (3.1) 

In this work we are interested in variations 5N which result from a change Sr(k) in 
the reflection amplitude. Due to the first condition in (2.2) we must impose 5r(— k) = 
6r*(k). Thus, with no loss of generality, we take the positive components r(k), k > 
as the independent functional variables. Keeping that in mind, we obtain from (2.7) 
and (2.9) that 

SN(a,b) _ 1 ik(a+b) 
5r(k) ~ 2n 6 ■ {6 - Z) 
Thus, from (3.1) and (3.2) we obtain 

5K(x,y) JSK(x,y) 6N(a,b) 



8r(k) J 8N(a,b) 5r{k) 

^e iky + J R(y,a;l)e lka dc}j [e ikx + J K(x, b) e ikb db^j . (3.3) 

From (2.4) and from the reality of K(x, y) and k, we recognize the last factor in (3.3) 
simply as ?p*(x, k). Thus, 



SK(x,y) _ 1 
Sr(k) 2tt 



e iky + 



J R(y,a;l)e ika da\ip*(x,k) (3.4) 



In order to simplify (3.4) further, we have to study the function 

oo 

uj(y,x;k) = J R(y,a;l)e ika da. (3.5) 

X 

Observe from (2.10) (considered with a generic given function f(y)), that fl(y, x; k) = 
e ik y _|_ u ^y^ x . ^ j g un iq Ue so l u tion of the Fredholm equation £l(y, x; k) = e tky + 

oo 

/ N(y, z)Q(z, x; k) dz, from which we infer that u(y, x; k) is the unique solution of 

X 

oo 

to(y, x; k) = —G(y, x; k) + J N(y, z)lo(z, x\ k) dz (3.6) 

X 

where 

G(y, x- k) = [°° F(y + z)e tkz dz . (3.7) 

Jx 

Thus, d x u>(y, x\ k) satisfies 

= F{x + v) ( e - + fc) ) + JjVfe, ,3 . (3.8) 

X 

From (3.5), (2.13) and the fact that R(y, z; 1) = R(z, y; 1), we obtain that u(x, x; k) = 

oo oo 

/ R(a, x; l)e lka da = J K(x, a)e lka da . Thus, the inhomogeneous term in (3.8) is sim- 

X X 

ply F(x + y) (e tkx + J K (x , a) e lka da^j = F(x + y)ip*(x,k) = — f(y)tp*(x, k), where we 
used (2.4). It is the given function f(y) in (2.8) multiplied by a ^-independent factor 
—ift*(x, k). Thus, from linearity, the unique solution of (3.8) is simply the solution of 
(2.8), multiplied by the same factor, namely, 

MV d *' k) =-K(x,y)r(x,k). (3.9) 

oo 

The initial condition for this equation at x = y is obviously co(y, y;k) — J K(y, a)e tka da 

y 

ip*(y,k) - e iky . Thus, 

y 

u{y, x;k)= f ^j*{z, k)K{z, y) dz + if)*{y, k) - e iky . (3.10) 



X 



Substituting this result into (3.4) we obtain our first main result: 

^^ = -i- L*(y,k)+fr(z,k)K(z,y)dzU*(x,k). (3.11) 



Since K(x,y) is a real kernel, we can write (3.11) alternatively as 

SK(x,y) 



Sr*(k) 



~ U(y,k) + Ji;(z,k)K(z,y)dzU(x,k). (3.12) 



The formula for S fJ?^) * s ^ e k e Y ^ or obtaining the functional derivatives of the wave 
function ip(x,k) and potential V(x) with respect to the reflection amplitude r(k), 
since the former are linear in K(x,y). Thus, from (2.5) and (3.12) we obtain our 
second main result: 

SV ( X ) = _ 2 — 5 J^^A = I A. tfU k ) . (3.13) 
5r*(k) dx 5r*(k) n dx 

Similarly, from (2.4) and (3.12) we obtain 

Sip(x,k) 1 



5r*(q) 2n 



e~ iky U(y,q) + J iP(z,q)K(z,y)dz\ dy 



i>(x,q). (3.14) 



Reversing the order of integrations in the second integral according to 

oo / y \oo /oo \ 

J e- lk y[j^z A )K{z,y)dz\ dy = J if>(z,q) I J K \z , y) dy) dz 

x \x I x \z / 

and recognizing the y-integral on the right hand side of the last equation as 

oo 

j e ~ lky K(z, y) dy = ip(z, k) — e~ lkz , we obtain from (3.14) our third main result: 



5ip(x, k) 



~{j i/i(z, q)iP(z, k) dz) V>(z, q) . (3.15) 



Sr*(q) 

It is interesting to note that both (3.13) and (3.15) are expressed purely in terms of 
the wave function ip(x,k). Note, however, that (3.13) is local in ■?/>, whereas (3.15) is 
highly nonlocal. It would be interesting to interpret these features from a physical 
point of view. 

4 Consistency Checks of (3.13) and (3.15) 

The integrated form of (3.13) should agree with the derivative Sr * (k ) J V(x) dx ob- 

— oo 

tained from the trace identity (1.5). Note that (1.5) is expressed purely in terms 



of the positive Fourier modes of r(k), the independent functional variables in our 
problem. Thus, taking the derivative of (1.5), we obtain 

2 r(k) 



5r*(k) 



J V(x) 



dx = — 



7r Sr*(k) 



log 1 - \r(q)\ dq 



*\t{k)\< 



(4.1) 



where we used \t(k)\ 2 = 1 — |r(A;)| 2 . This result should be confronted with the inte- 
grated form of (3.13). We obtain from the latter 

5 



1 lim 



, V(x) dx — - 



ip (x, k) 



(4.2) 



From (1.8) we observe that ip(L,k) ~ e %kL . Similarly, from (1.10) and (1.4) we 
deduce the asymptotic behavior 

1 



4>(-L,k) 



£ ikL _ r (k) c -ikL 



t*(k) t(k) 

Substituting (4.3) and tp{L,k) ~ e~ ikL in (4.2), we obtain 



(4.3) 



Sr*(k) 



°r 1 

/ V{x)dx = - hm 

J vr r,_.. 



L— >oo 



e~ 2lkL 1 



t 2 



_J_ 2ikL ^_ 

(t*y \t\ 2 



(4.4) 



which coincides in the limit with (4.1) in the sense of distributions, since the first two 
rapidly oscillating terms on the right hand side of (4.4), when smeared against any 
continuous bounded test function U(k) will integrate to zero in the limit L — > oo, due 
to the Riemann-Lebesgue lemma. Thus, (3.13) has passed its first test. 

A less trivial test of (3.13) arises from comparing the complex conjugate equation 
of (3.13), = ip* 2 (x, k) and (1.2). We see that one kernel is the inverse of the 
other: / ™ 



8V(x) 



dx = 5(k — q) . Thus, we must verify that 



d 



(4.5) 



is equal to 8(k — q). Using the identity 



, d 

dx 



G dx 



CO 

(FG) 2 ] _ oo + / FG {Fd x G - Gd x F) dx 



with F = (f)(x, k) and G = ip*(x, q), we write (4.5) as 
t 2 (k) (l 



oo 



T ( k ^) = ^i l-[(<P(x,k)r(x,q)f}™ oD + I ( />(x,k)r(x,k)W(x,k,q)dx\ (4.6) 



-oo 



where 

W(x, k, q) = <f>(x, k)d x ip*(x, q) - i^*(x, q)d x <f)(x, k) (4.7) 

is the Wronskian of 4>(x,k) and ip*(x,q). Next, from the Schrodinger equation for 
these two functions it is easy to obtain the relation 

4>{x, k)r(x, q) = —^d x W(x, k, q) . (4.8) 
k z — q z 

Substituting the last equation in (4.6) we see that the integral is given entirely by 
boundary terms as 

m 



T(k,q) 



(<f>(x, k)r(x, q)Y + T2 ~ 2 W 2 (x, k, q) 

k z — q z 



(4.9) 



<iiiik 

In order to proceed in the clearest and simplest way we shall make a few observations. 
First, observe from (1.4), (1.8) and (4.3) that the difference (<j>(L, k)ip*(L, q)) 2 - 
(0(— L, k)vf)*(—L, q)) 2 is rapidly oscillating in the limit L — > oo for all k, q. Thus, as a 
distribution acting on bounded continuous functions of k and q it tends to zero. Thus, 
we should focus on the second term. When k 2 — q 2 ^ 0, the difference W 2 (L, k, q) — 
W 2 (—L, k, q) oscillates rapidly, and tends to zero, in the sense of distributions, in the 
limit L — > oo, similarly to the first term. Recall that for k 2 = q 2 , i.e., when <p(x,k) 
and ip*(x, q) correspond to the same energy, W(x, k, q) is a constant, and thus cannot 
oscillate. In that case, however, the denominator in front of this term vanishes. Thus 
we should study the limit q 2 — > k 2 carefully. 

Let us concentrate then on the region k 2 ~ q 2 . Since we have taken the positive 
Fourier modes of r(k) as the independent functional variables, it is enough to study 
the case q — > k. Thus, assuming k — q = e with ||| << 1 and setting k + q ~ 2k, we 
obtain from (1.4), (1.8) and (4.3) that 



W(L,k,k-e) ~ -i- (2ke~ ieL - er(k)e 



2ikL 

t(k) 



W(-L, k,k-e) ~ \2ke* L - rf(fc)^e^J . (4.10) 
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Substituting (4.10) and k 2 — q 2 ~ 2ke in (4.9), and recalling the first observation 
made right below (4.9), we finally obtain 

r(t,,) = fE^ + H.O.T., (4.11) 

where the acronym R.O.T. stands for the rapidly oscillating terms that tend to zero 
in the sense of distributions. The first term, tends of course to 5(k — q) as L — > oo, 
as desired. (3.13) has passed the second check! 

Our last consistency check concerns (3.15). The point is that for x — > — oo, the 
integral on the right hand side of (3.15) is the orthogonality relation 

4>(z, q)4{z, k) dz = (5(k + q)- r(k)S(k - q)) . (4.12) 

Let us sketch the proof of (4.12): in analogy with (4.8) we deduce that 

V>(x, k)ifj{x, q) = —L^d x W{x, k, q) (4.13) 
k — q A 

where 

W(x, k, q) = ifj(x, k)d x ifj(x, q) - ip(x, q)d x ip(x, k) (4.14) 
is the Wronskian of the functions involved. Then, after integration, we obtain 

J 4>(z, qMz, k) dz = ^—^ [W(L, k, q) - W(-L, k, q) 

and consider this result as L — > oo, using the asymptotic behavior (1.8) and (4.3). For 
generic values of k, q, the difference of Wronskians is a rapidly oscillating function, 
which as a distribution, tends to zero as L — > oo. As in the previous discussion, we 
observe that when k 2 = q 2 , W(x, k, q) is independent of x, and thus does not oscillate. 
By studying the limit q — > ±k carefully, we deduce (4.12). 
On the left hand side of (3.15) we have 

5r(-q) 5r(-q) [t(-k) e t{kf ) ' ^ 



11 



where we used (4.3). Then, from (2.3) it follows that 




2nit(k) \t(p) | 2 \p-k-ie p + k + ie) 



( 



p — k p + k 



1 1 



) 



(4.16) 



Applying (4.16) to (4.15), and studying the resulting expression around k ~ q > 0, 
we find (keeping only the singular terms) ^ ^^ 2 ijj(—L, k), in accordance with the 
right hand side of (3.15) and (4.12). 

Appendix: Miscellaneous Technical Details 

A.l Perturbation Theory in the Continuum: Proof of (1.2) 

Consider a perturbation 8V(x) of the Schrodinger equation (1.3). To make the prob- 
lem well-defined, we assume that SV(x) is localized around some point x . 

Under this perturbation, the solution of (1.3) will be shifted ty(x) — > ^(x)+8^(x), 
and our task is to compute 8^f(x) to first order in 8V(x). To this order, we have to 
solve the equation 



The general solution is given in terms of the Green's function G(x, y\ k) of the operator 
H — k 2 (defined with the appropriate boundary conditions) as 



Since this expression vanishes for 8V(y) = 0, there is no term which is a solution of 
the homogeneous equation. 

The proper Green's function is the one which decays exponentially whenever one 
of its coordinate arguments tends to ±oo, when k is lifted to the upper half complex 
plane. There is a unique such Green's function: 



(H - k 2 )8^(x, k) = -8V(x)V(x, k) . 



(A.l) 




(A.2) 



G(x,y;k) 



1 



[6 (x - y) ip*(x, k)(f)(y, k) + 9 (y - x) tp*(y, k)<f>(x, k)] 



(A.3) 



w(k) 



12 



where 



w (k) = -W(x,k,k) 



2k 
it(k) 



(A.4) 



from (4.7). Thus, in particular, for ty(x, k) = (f>(x, k) we obtain 
1 



8<f>(x, k) = ~ i(;*(x, k) J (j)(y, k) 6V(y)<f>(y, k) dy + <j>(x, k) J i(j*(y, k) SV(y)4>(y, k) dy 

— OO X 

(A.5) 

In the limit x — > — oo, we see that 

(p(x,k) 



5(f>(x, k) 



w 



Thus, the function 



4>(x, k) 



^*(y, k) 8V(y)<f)(y, k) dy . 



(x, k) + 5(f>(x, k) 



(A.6) 



1-i / r{y,k)8V{y)4>{y,k)dy 



(A.7) 



tends to e~ tkx as x — > — oo, and should be identified with the "0-function" (1.4) of 
the perturbed potential V(x) +SV(x). In (A.7) we should keep, of course, only terms 
up to linear order in 5V(x). As x — > oo, we see from (A.5), (1.4) and (1.8) that 
4>(x, k) — > j^ye~* fe:r + j^e lkx , where t(k) and f(k) are the new scattering amplitudes, 
given by 



1 


1 


t(k) 


t(k) 


r{k) 


r(k) 


t(k) 


t(k) 



1 + 



1 00 

- I r(y,k)5V(y)<j ) (y,k)dy 



-oo 



(A.8) 



from which we infer, after some simple algebra, involving (1.9) that 

oo 

<5r(fc) = f (fc) - r(k) = — J (t(k)<P(y, k) f 8V(y) dy , (A.9) 

— oo 

from which (1.2) follows. 
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A. 2 Reflect ionless Potentials 

We demonstrate in the following the positivity of the Fredholm determinant of the 
GLM equation (2.6) in the case of reflectionless potentials. 

For reflectionless potentials, where r(k) = 0, N(y, z) — — Ya=\ cfe~ Kl( - y+z ^ in (2.8) 
is a degenerate kernel of finite rank N. Thus, the integral equation degenerates into a 
system of linear equations, and the Fredholm determinant becomes an N dimensional 
determinant. For this reason, reflectionless potentials are so easy to treat within the 
formalism of 1ST. 

The Fredholm determinant (at A = 1) is given by (ci • • • c^) 2 det A, where 

A mn = 5 mn + -E^H- e -i^> . (A. 10) 

This matrix is a member of a family of matrices of the general form 

Aiy) — x I VmVn (AlU 

where v n are the components of a real vector and v a real number. (Recall also that 
all K n > 0.) The matrix (A. 10) corresponds, of course to v n = c n e~ KnX and v — 1. 
It is easy to demonstrate that (A. 11) is a positive definite matrix for any v > 0. To 
prove this, it is enough to verify that £ T A£ > for any real vector £. Thus, consider 

t, At; - ^ ^ + 2^ , y 
n,m K^n ~r rv m ) 

1 f oo ( N \ 2 

= Z T t+—j o I E^^e" K " s 1 ds>0 (A.12) 

manifestly. Thus, in particular, (A. 10) is positive definite and so is the corresponding 
Fredholm determinant. 
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